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We prove combinatorially the explicit relation between genus filtrated s-loop means of the 
Gaussian matrix model and terms of the genus expansion of the Kontsevich-Penner matrix model 
(KPMM). The latter is the generating function for volumes of discretized (open) moduli spaces 
Mg's C given by N 9tS (Pi,... ,P S ) for (Pi,... ,P S ) G Z+. This generating function therefore enjoys 
the topological recursion, and we prove that it is simultaneously the generating function for an¬ 
cestor invariants of a cohomological field theory thus enjoying the Givental decomposition. We 
use another Givental-type decomposition obtained for this model by the second authors in 1995 in 
terms of special times related to the discretisation of moduli spaces thus representing its asymptotic 
expansion terms (and therefore those of the Gaussian means) as finite sums over graphs weighted 
by lower-order monomials in times thus giving another proof of (quasi)polynomiality of the dis¬ 
crete volumes. As an application, we find the coefficients in the first subleading order for in 

two ways: by using the refined Harer-Zagier recursion and by exploiting the above Givental-type 
transformation. We put forward the conjecture that the above graph expansions can be used for 
probing the reduction structure of the Delgne-Mumford compactification A4 giS of moduli spaces of 
punctured Riemann surfaces. 


1 Introduction 

After Harer and Zagier derived m the celebrated recursion formula describing means of traces (tr H k ) 
of powers of the N x A-Hermitian matrix H over the Gaussian ensemble, progress in this direction was 
mainly due to explicit (commonly multiple-integral) formulas for the generation functions for multi¬ 
trace means of the form (H§L» , tr H ki ) conn . Brezin and Hikami m used the replica method ameliorated 
in |35| for obtaining exact s-fold integral representation for these quantities valid for all N. There has 
been revived interest in the multi-trace means due to their relation to topological recursion [10L [20lj and 
quantum curves mmm One aim of the present paper is to describe consequences of these recent 
ideas. In particular, the relation to the Euler characteristic of the moduli space of curves, which is 
the main application in [27] , is shown here to go much deeper and culminates in a cohomological field 
theory [33] equivalent to the multi-trace Gaussian means. 

In this paper we exploit the interplay between three main tools for studying the multi-trace Gaus¬ 
sian means: 

1. A matrix model, known as the Kontsevich-Penner matrix model (KPMM), with external field 
that describes all multi-trace Gaussian means, introduced in J2j and developed further in |3]. 
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2. Exact formulae, or equivalently integrals over N x N Hermitian matrices for fixed N, in place 
of asymptotic expansions of means in N which gives rise to surfaces and their genus. 

3. Topological recursion satisfied by resolvents storing the multi-trace Gaussian means, which leads 
to a cohomological field theory with primary correlators given by x(Af SjS ). 

The construction of the Kontsevich-Penner matrix model uses a standard combinatorial approach. 
As was shown in m using the Virasoro constraints and in [33j using the direct determinant relations, 
this model is equivalent to the Hermitian matrix model with the potential whose times (coupling 
constants) are related to the external-matrix eigenvalues via the Miwa-type transformation and whose 
matrix size is the coefficient of the logarithmic term. More precisely, the statement, which appears 
here as Corollary [U is that the free energy for the KPMM is a primitive (antiderivative) for the 
resolvents of the Gaussian matrix model. The resolvents storing the multi-trace Gaussian means are 
naturally described as meromorphic (multi)differentials with zero residues over a rational Riemann 
surface, known as the spectral curve , hence their primitives are meromorphic functions on the spectral 
curve. These primitives are conjecturally related (and proven in the Gaussian case [37] ) to the so 
called quantum curve which is a linear differential equation that is a non-commutative quantisation 
of the spectral curve. Corollary [T] can be used to prove the relationship between the spectral and 
quantum curves—it shows that the wave function arising out of the spectral curve is a specialization 
of the free energy for the KPMM which satisfies the second order differential equation that is the 
quantum curve. 

The geometric content of the Kontsevich-Penner matrix model also proved to be rather rich: its 
free energy was related to structures of discretized moduli spaces in [8]; more recently in [39] and [36] 
it was identified with the generating function for discrete volumes 1V S s (Pi,..., P s ) —quasi-polynomials 
introduced in [38] that count integer points in the interiors M g , s of moduli spaces of Riemann surfaces 
of genus g with s > 0 holes with the fixed perimeters Pj € Z+, j = 1,..., s of holes in the Strebel 
uniformization. Moreover, it was shown in [9] that, when being expressed in the times T± n originated 
from the discretizations of the moduli spaces, this model admits a decomposition into two Kontsevich 
models related by a Bogolyubov canonical transformation—to the best of our knowledge it was the 
first example of the Givental-type decomposition formulas [23] . 

We use the decomposition formula of [9] to express the terms P g%s of the free energy expansion 
of this model in the form of finite sums over graphs whose vertices are terms of the expansion of the 
free energy of Kontsevich matrix model, internal edges correspond to quadratic terms in the canonical 
transformation operator, external half edges (dilaton leaves) correspond to the constant shifts of the 
higher times, and external legs (ordinary leaves) carry the times T^ n . This graph representation 
provides another proof of quasi-polynomiality of N 9jS (Pi, ..., P s ) and that these quasi-polynomials 
depend only on even powers of Pi. 

From [19] and m we know that the terms of topological recursion EBMnamin based on 
a certain spectral curve satisfying a compatibility condition (relating the wo,i and wo ,2 invariants) 
describe ancestor invariants of a cohomological field theory, or equivalently a Frobenius manifold. 

A fundamental family of Frobenius manifolds described by Dubrovin are Hurwitz spaces. For 
g = (g±,... , g n ) define the Hurwitz space H ggj to consists of genus g branched covers of the sphere 
with n labeled points over oo of ramification profile (g \,..., g n ) and simple ramification over P 1 — oo. 
It has dimension \g\ + n + 2g — 2 where \g\ = g\ + • • • + g n . 

The 2-dimensional Hurwitz-Frobenius manifold i7 0 ,( 1 , 1 ) consists of double branched covers of the 
sphere, with two branch points and unramified at infinity. Its free energy is 

F = 7^*0,0*0,1 + ^o,l lo g*0,l 


(1.1) 
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Note that expression m appears as a standard term (the perturbative part) in the expansion of any 
matrix model upon identification of *o,i with the normalized number of eigenvalues. 

We discuss the combinatorial relation between the discrete volumes and Gaussian means using 
the cohomological field theory (CohFT) description, which relates the discrete volumes to ancestor 
invariants of a CohFT. These ancestor invariants are evaluated already in terms of the closed moduli 
spaces Mg yS compactified by Deligne and Mumford. We thus provide explicit description of these 
invariants in terms of the Gaussian means Ws 9 \x i, ..., x s ). 

In 0, M the conjecture was put forward that the very same KPMM expressed in times T 2n de¬ 
scribes the Deligne-Mumford stratification of the closed moduli spaces M.g )S while reduction coefficients 
Cg,s,r q [H] related to ancestor invariants in the above scheme were conjectured to be positive rational 
numbers proportional modulo some presribed combinatorial factors to the numbers of screens m 
for the corresponding Riemann surfaces. Although the original conjecture of [8], which assumed the 
existence of covering tori for M. g , S i turned out to be incorrect, this interpretation nevertheless helps 
to explain why the discrete volumes N g s (Pi, ..., P s ) are Z2 quasi-polynomials modulo the parity of 
Pi, not polynomials, and why the dependence is only on even powers of Pi'. In formulas for the moduli 
space multi-component reductions the two sets of times undergo synmretrization depending on the 
multicomponent reduction structure, and we therefore lack the single symmetrized expression when 
evaluating the inverse Laplace transform. 


The paper is organized as follows. In Sec. [2] we establish the equivalence between the Gaussian 
means (the correlation functions) and the terms of expansion of the KPMM free energy. 

We devote Sec. [3] to describing results of [8] and [9] concerning the open discrete moduli spaces by 
which we produce the formula relating the above Gaussian means and the discrete volumes in a purely 
combinatorial way. We also demonstrate how the quantum curve can be obtained as a specialization of 
the KPMM to the case of unit size matrices. We describe the Givental-type decomposition formulas 
for the KPMM obtained in [9] representing them in terms of graph expansions for the free energy 
terms. We use this graph representation for proving the quasi-polynomiality of the discrete volumes 
and for making a link to CohFT. Finally, we discuss a link provided by this representation to a 
Deligne-Mumford-type stratification of moduli spaces. 


In Sec. 01 we identify the Gaussian means expansion terms with the ancestor invariants of a 
cohomological field theory using the results of |16] and m- The decomposition thus obtained has 
a canonical Givental form differing therefore with the decomposition in the preceding section. We 
show that the coefficients of this (“alternative”) decomposition, or equivalently, the coefficients of the 
Laplace transforms of the quasi-polynomials N gyS (P\,... ,P S ), are the special coefficients bj?\ which 
are also the coefficients of P g (k) from [5l IB], that represent in the “most economic” way to record 
the genus filtered s-loop means Ws 9 \x 1 ,... ,x s )—for example, for s = 1, we need only g coefficients 
b^ 9 \ k = 0,... g — 1. We identify the coefficients \) 9 ' 1 with simple linear combinations of the CohFT 
ancestor invariants. 


In Sec. [H we concentrate on the case of a one-loop mean satisfying the Harer-Zagier (HZ) recurrence 
relation which is proven using integrals over N x N Hermitian matrices for fixed N. On the basis 
of the original HZ relation, we obtain a more effective recursion relation for the coefficients l)f ', 
k = 0,..., g — 1, of shapes and formulate the recursion procedure for obtaining the higher coefficients 
b^ 9 l s for a fixed s at all g. We prove that b^ 9 \ k = 0,..., g — 1, are positive integers for all g and 
k, thus providing further evidence for the main conjecture of [5], which would imply this positivity 
property. We present two alternative derivations of the first sub leading coefficient b[ 9 \ 2 : from the HZ 




recurrence relations and from the decomposition formulas of Secf3l 


2 The matrix-model representation for the multi-loop Gaussian means 

We consider a sum of connected diagrams with s backbones, or loop insertions, each carrying the 
corresponding variable Ui, i = 1,..., s. Our aim is to formulate the matrix model that describes all 
genus-g contributions in terms of shapes —the connected fatgraphs of genus g with s faces and with 
vertices of arbitrary order greater or equal three; according to the formula for the Euler characteristic, 
for a fixed g and s we then have only a finite number of such graphs, and we let denote this 
finite set. In fact, r g>s enumerates the collection of cells in the canonical (Strebel-Penner) ideal cell 
decomposition of moduli space -M g , s , the combinatorial part of which we shall recall here. 

The origin of the word shapes here come from [j5j, where we considered the Poincare dual graphs 
of such fatgraphs, and by analogy with the arguments in that paper, one can easily see that is in 
bijection with circular chord diagrams which are also ’’shapes” in the terminology of [5], that is chord 
diagrams which are seeds and which has no one-chords, again in the terminology of [5|. 


2.1 Multi-loop Gaussian means 

We begin by considering the correlation functions, or means, 
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where T~Ln is the set of Hermitian N x N matrices. By Wick’s theorem, any correlation function 
of the form (12.21) can be calculated as the sum over all possible (complete) pairings between matrix 
entries My, and every pairing is given by just the two-point correlation function (HijHkj) = jj5u5jk- 
Such pairings gives rise to fatgraphs containing ordered set of s vertices of valencies ki, i = 1,..., s, 
and Yll=\ ki/2 edges corresponding to the pairings (represented by double lines, thus inducing the fat 
graph structure). Furthermore for each vertex we also have a first incident edge given. We denote this 
set of fatgraphs r(fci, ..., k s ). Then the above sum reads 
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where 6 ( 7 ) is the number of boundary components of 7 . 

Let r(& 7 ,..., k s ) c be the subset of r(& 7 ,..., k s ) which consist of connected such fatgraphs and we 
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Second, we consider the 1/iV-expansion of the connected correlation function 
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to segregate its part corresponding to connected fat graphs of genus g. Wick’s theorem therefore 
implies that 

/T-r i \ conn ^ 
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where T g (k\, ..., k s ) c is the subset of r(&i,..., k s ) c consisting of genus g fatgraphs. We now consider 
the union 

F C g, S = U f 

{ki,...,k a }ez s + 


By blowing up each vertex to a backbone we exactly get the genus g chord diagram on s linear 
backbones as considered in 0 , 0 (see Fig. HD- 




Figure 1: Transition from ciliated multivalent vertices to linear backbones. 
We have the following formula 
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where T g (k \,..., k s ) c is the set of connected fat graphs of genus g with s (nonciliated) ordered vertices 
of valencies k\,... ,k s and Aut ( 7 ) is the automorphism groups of the fatgraph 7 which respects the 
ordering of the vertices. We obtain the actual multi-loop case (up to overall sign) by differentiation: 
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By combining formula (12.51) with (12.71) . we find that 
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In order to proceed to the Poincare dual fatgraphs, which provides an equivalent description of 
chord diagrams in terms of shapes understood as closed fatgraphs with vertices of order not lower than 


















three, it is more instructive to consider the non-ciliated fat graphs and perform the proper resummation 
(see the next subsection). We observe that the effect of differentiating w.r.t. the variables Ui as in 
(12.71) is precisely to add a cilium at vertex i. 


2.2 Summing up rainbow and ladder diagrams—formulating the matrix model 

In the expression (12.61) . we may perform a resummation over planar diagrams. To this end, we shall 
refer to a planar chord diagram on an interval as a rainbow diagram as depicted in Fig. [2j The number 
of rainbow diagrams of course have the generating function 
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so summing over rainbow diagrams is thus mere Catalan number counting, which effectively reduces 
to replacing the original edge of a chord diagram by a thickened edge carrying the factor f{uf). 
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Figure 2: The procedure of summing up rainbow diagrams of chords (dashed lines) for a single backbone (solid 
lines). The result is the new edge of the backbone indicated by a double line (the thickened edge). 

We now consider the summation of ladder-type diagrams, where a “rung” of the ladder joins two 
cycles that carry in general distinct (but maybe coincident) indices i and j (see an example in Fig. [3]). 
Each ladder contains at least one rung, which is a chord carrying the factor UjU-j. We obtain the set of 
new edges and new vertices if we further blow up cycles of double-line backbone edges until they will 
be joined pairwise along rungs (each containing at least one rung); disjoint parts of these cycles will 
then constitute loops of lengths 2r^ > 6 alternatively bounded by r & rungs (the chords) and rf. double 
lines (the thickened edges of circular backbones, possibly with different backbone indices of different 
edges); these loops then become vertices of the respective orders > 3 of the new fat graph. 



Figure 3: The procedure of performing a sum over ladder diagrams. The thickened edges associated with the 
selected ladder are painted dark. The crosshatched domain will become a four-valent vertex in the new fat 
graph, the white domain to the left of the crosshatched one becomes a three-valent vertex, etc. Ladders turn 
into edges of the new fat graph. The double solid line in the right-hand side of the figure is the edge obtained 
from darkened thickened edges in the left-hand side of the figure. 

The diagrammatic expression (Fig. [3]) demonstrates that we have exactly one factor f{ui) and 
one factor f(uj) per each rung of a ladder (we introduce the same cyclic orientations on all circular 
backbones): these are the factors associated with those thickened edges of backbones that precede the 















points of attachments of the rung to the ith and jth backbones. This assignment is thus uniquely 
determined and no thickened propagators remain unassigned, so the corresponding sum of ladder 
diagrams reads 
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We arrive at the main result of this section: 
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Theorem 1. The genus-g term of the (nonciliated) s-backbone case is given by the following 
(finite!) sum over fatgraph shapes T ffiS of genus g with s edges whose vertices have valences at least 
three: 
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:= Fg 9 \X ±,..., X s ), ( 2 . 11 ) 


where ± denotes the two sides (faces) of the edge e. The quantity Fg 9 \ Ai,..., X s ) in the right-hand 
side is the term in the diagrammatic expansion of the Kontsevich-Penner matrix model j!2jf given by 
the following normalised integral over Hermitian N x N-matrices X: 
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where the sum ranges over all stable curves (2 g + s > 2) and A is a diagonal matrix with the numbers 
e Xi on the diagonal. 


The proof beyond the earlier works [12] . [42] has essentially been given already. Let us further 
simply observe that in replacing boundary cycles in ladder diagrams by punctures and rungs by arcs 
connecting those punctures, we reproduce the small elaboration of the original arc families in m where 
parallel copies of arcs are permitted in the current instance. The underlying idea is then completely 
obvious graphically, where (possibly empty) copies of boundary-parallel rainbow diagrams are here 
permitted in between arcs. Notice that the old dual fatgraph to the arc family in the original sense of 
[42] is the new fatgraph of the held theory defined above. 

Differentiating the relation p.llft w.r.t. A i on the right-hand side we obtain the standard loop 
means, or (connected) correlation functions Wg 9 \x i,... ,x s ), x % = e Xi + e~ Xi , of the Gaussian matrix 
model enjoying the standard recursion relations of the topological expansion |T8] . m- As an immediate 
consequence of the previous result, we have 

We introduce the standard definition of the connected correlation functions (the s-loop means, or 
resolvents) of the Gaussian matrix model with the potential ^-tr H 2 : 
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Although the difference between thus defined quantities and original means (12.61) as well as with the 
generating functions C g ^ s (z) for s-backbone chord diagrams of [5j and [6] is only in changing the 
normalisation (see Remark [3]) precisely these resolvents (12.131) turn into symmetrical s-forms in the 
formalism of topological recursion of [18] . [10] thus playing an instrumental role in this formalism. 


Corollary 1. The exact relation between resolvents 12.131) and the terms of the free-energy ex¬ 
pansion of the Kontsevich-Penner matrix model reads 
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The quantities Ws(x i,..., x s ) here enjoy the standard topological recursion II Of . J7f for the spectral 
curve x = e A + e _A , y = ^(e A — e~ A ). 


3 Matrix models and geometry of moduli spaces 

3.1 The Kontsevich construction for evaluating intersection indices 


We begin with recalling the cell decomposition of moduli spaces of Riemann surfaces of genus g with 
s > 0 marked points which was proved independently by Harer [26], who presented a proof of Mumford 
using Strebel differentials [J3]> an< 4 by Penner m using hyperbolic geometry. This cell decomposition 
theorem states that strata in the cell decomposition of the direct product M g , s x R+ of the open 
moduli space and the s-dimensional space of strictly positive perimeters of holes are in one-to-one 
correspondence with fat graphs of genus g with s faces (those are the shapes from Sec. [2]) whose edges 
are decorated with strictly positive numbers U £ R+. The perimeters Pj, I = 1,..., s are the sums of 
li taken (with multiplicities) over edges incident to the corresponding face (boundary component, or 
hole). So it is natural to call them the lengths of the corresponding edges. 

One associates the Chern classes to the corresponding holes, and Kontsevich m found a remark¬ 
ably simple formula for representatives of these Chern classes i/jj: 
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where Pi = ^2j=i hj an d h,---,ini enumerates edges incident to the I’th boundary component, 
which is therefore an ni -gone and these edges are linearly ordered in accordance with the orientation. 
Changing the starting edge results in adding an exact 2-form to (13.151) thus not affecting the result of 
integration. The quantities to be calculate are the intersection indices 
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which do not depend on actual values of Pj being purely cohomological objects. The useful tools that 
enables the determination of the generating function for these numbers are, first, to multiply every ij> < j I 
by Pj dl and, second, to perform the Laplace transformation w.r.t. all Pj, that is, to integrate over the 
total direct product A 4 gyS x M^_. All unwanted factors arising when the external derivative d acts on Pi 
in the denominators are then canceled by the integration measure dP\ A • • • A dP s , and only the total 
volume form coming with the factor 2 2ff ~ 2+s remains in the integral, which therefore transforms into 
the diagrammatic expansion of the Kontsevich matrix model; this procedure was extensively reviewed 
















in the literature (see, e.g., [32]). On the other hand, this integration merely gives 
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The left-hand side of (I3.16P can be presented as the sum over three-valent fat graphs with the 
weights 1/(A/ 1 + A/ 2 ) on edges where I\ and 1-2 are indices of two (possibly coinciding) cycles incident 
to a given edge. Also a factor 2I appears (where | V| and |L| are the cardinalities of the respective 
sets of vertices and edges). The generating function is then the celebrated Kontsevich matrix model 
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are the times of the Kontsevich matrix model. 
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3.2 The discretization of moduli spaces 

The material in this subsection (and only in this subsection) is not rigorous and we need it only 
as a motivation for what to follow. We suppose here the existence of a discretised version of the 
Kontsevich construction in Subsection 13. II and that such a hypothetical theory would also make sense 
on the boundary. 

As was proposed in [8], set all the lengths of edges of the Penner-Strebel graphs to be nonnegative 
integers k € Z + , i = 1,..., |L| < 6g — 6 + 3s. 

Were a proper cohomological theory of discretely acting operators to exist, we could introduce 
special representatives of the corresponding Chern classes Vfr i n the spirit of the original Kontsevich 
construction m introducing also the analogue of the 12 two-form 


n = Y,Pi$i- (3-19) 
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We would then replace integrations w.r.t. Z, by half-infinite sums over integer points; in [8] we 
called these lattices endowed with discrete integration the discrete moduli spaces 

When integrating the above forms over the discrete moduli spaces, as in the case of the standard 
moduli spaces, the integral /^disc 11/ Pj ‘'^i' does not make sense until we (i) include the closure of the 
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moduli space and (ii) perform the discrete Laplace transformation w.r.t. the perimeters Pj weighted 
with the standard function Pl> ' 1 thus reconstructing, as in the original Kontsevich approach, the 

total volume form on the totally discretized space x Z^_. Then, as in the Kontsevich case, we 

have the purely combinatorial relation between the highest-order forms, 
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which does not depend on the numbers di (provided dj = 3g — 3 + s) and on the combinatorial 
type of the fat graph on whose edges have length li. 







An important subtlety is that since the length of every edge enters twice in the sum ^j =1 Pj, this 
sum is always a positive even number, and we must take this restriction into account when performing 
the discrete Laplace transformation. 

Because the Deligne-Mumford compactification of the moduli space is purely combinatorial, prop¬ 
erly defined integrals over continuous and discrete moduli spaces must match each other, and instead 
of the formula (13.1611 we presumably obtain a formula of the type 
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where the new times 
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are discrete Laplace transforms. 

In the next subsection, we demonstrate that the KPMM (12.1211 corresponds to integrations over 
the open discrete moduli spaces. Note here that following the Deligne and Mumford approach, we can 
represent integration over the open moduli spaces M g>s (discrete or continuous) using the stratification 
of closed moduli spaces, as an alternating sum over reduction components, 
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where |r| is the reduction degree equal to 3g — 3 + s — Yhftgi — 3 + si) and c%’ s are positive rational 
numbers counting, roughly speaking, the numbers of the given type of reduction (boundary strata) 
per copy of M. g , s in the corresponding cell decompsition of M g , s (also related to the number of screens 

Would the Kontsevich symplectic structure have an analytic continuation to the boundary, we 
might expect that we can use (|,3.22ll in order to present an “integration” of the form (see (|3.1911 1 
over an open moduli space M. g , s in the form of alternative sums of products of integrations over closed 
moduli spaces thus obtaining that 
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(3.23) 


where (r^, • • • Td,) are the standard Kontsevich intersection indices and we segregate ^-classes on 
reduction components into two groups: the original ^-classes (their total number is always s) and those 
obtained via reductions (we put dilaton leaves D^ dk in correspondence with these classes; originally 
all those classes were assumed to be just the trivial classes tq). Note that the symmetrisation w.r.t. 

















T + and T occurs in every reduction component separately, so the result contains mixed terms of T + 
and T~ thus being a quasi-polynomial. 

Although the naive form of (13.231) does not hold, notably, a factorisation of T± into expressions 
independently symmetrized in each reduction component fails and we have to introduce relations 
between dilaton leaves at different components of the reduction, thus producing (inner) edges of 
a special graph representation in Subsec. 13.41 the genuine resulting expression in terms of graphs 
presented in this subsection manifests at least some of the properties we might expect from (13.231) : 
it depends only on even times T^(A) and it exhibits an alternating sign structure w.r.t. the total 
powers of monomials in these times (see Remark 0]), so we expect that this model contains some useful 
information about the cell decomposition of the closed moduli spaces, and not only that of the open 
moduli spaces, on which we concentrate in what follows. 


3.3 Open discrete moduli spaces and the Kontsevich—Penner matrix model 


No matrix-model description of integrations over the closure of the moduli spaces exists (and would 
presumably just be obtain from the original Kontsevich intersection indices weighted by different 
times). On the other hand, integrals over discrete open moduli spaces are well defined for volumes of 
discrete open moduli spaces 

Following [38] we thus define the discrete volumes N gtS (Pi, ■ ■ ., P s ) which is a weighted count 
of the integer points inside x for fixed positive integers Pj, I = 1, ...,s, which are the 

perimeters of the holes (cycles). These numbers are equal (modulo the standard factors of volumes 
of automorphism groups) to the numbers of all fat graphs with vertices of valencies three and higher 
and with positive integer lengths of edges subject to the restriction that the lengths of all cycles are 
fixed. Using the above identity 
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O )) 


i=i 




where l e is the length of the eth edge and and are the indices of two (possibly coinciding) 
cycles incident to the eth edge, we obtain that 
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(3.24) 


We recognize in ([3.2411 the genus expansion of the KPMM (12.121) . We thus immediately come to the 
following statement: 


Lemma 1. The generating function for the Laplace transformed discrete volumes N gs (P\,,P S ) 
is the KPMM (12.121) . The correspondence (13.241) is given hy the formula 


v' n 2 ~ 

e K g ,s,p j ez + JV 


2 9a 2 - 2 9 - s i Vfli ,(P 1) ...,P.)e— £J= 1 p Pi _ f DXe~ aNtr [|A*A*+log(l-X)+x] 


fDXe- aNtr ^ AXAX -l x ^ 

where the sum ranges over stable curves with 2g — 2 + s > 0 and only positive perimeters Pp 


(3.25) 


Remark 1. The formula ([3.251) is valid at all values of N and A;. Specializing it to the case N = 1 
(when we have just an ordinary integral instead of the matrix one) and setting A; = A, a = 1 /h, and 
x = e A + e~ A , we obtain 


0 K,s, Pi ez + h 29 +‘- 2 N g KP?,.,P*)e~ X ^i= i p i _ 
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where the function 


OO 


F{h, x ) : = 

satisfies the second-order differential equation 


dt e ~\d 2 /2+ xt +\°gt) 


. ,9 


tr-zr-n + xh— + (1 — K) F(h, x ) = 0. 
ox A ox J 


We thus reproduce the equation of the quantum curve from 


Remark 2. A formal limit Pj —>• 0 in the Laplace transformed expressions correspond to A/ —>• 

—oo in (I3.25|) : in this limit, we obtain that A := diag({e A/ }) —>• 0, and 1V 9 )S (0,..., 0) becomes the 
(< 7 , s)-expansion term of the Penner matrix model [41] 

J DX e ~ aNtr [ lo g( 1 “ x )+ j,f ] 
fDXe - aNtrX2 / 2 

counting virtual Euler characteristics of moduli spaces (at negative a), i.e., 1V 9 iS ( 0,..., 0) = (— l) s x(-Mg,s)- 


Note that the discrete volumes are quasi-polynomials: their coefficients depend on the mutual 
parities of the P/’s and we present one more proof of this fact below (see Corollary [2]) . Because the 
corresponding generating function (12.121) is in turn (see formula (|2.14p ) related to the standard s-loop 
Gaussian means, or correlation functions wi 9 \ we have the following lemma. 

Lemma 2. The correlation functions wi 9 \x i,... ,x s ) of the Gaussian matrix model subject to 
the standard topological recursion based on the spectral curve x = e A + e _A , y = ^(e A — e _A ) are 
related to the discrete volumes by the following explicit relation: 


W^ 9) (e Al + e" Al ,..., e As + e" As ) = J] 
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N gt8 (P 1 ,...,P 9 ). (3.26) 


Example 1. We begin by considering A r 9 ) i(P), which are polynomials of degree 3g — 2 in P 2 (half 
the real dimension of the highest dimensional cells in the moduli space M g . i), are nonzero only for 
even P, and must vanish for all P = 2,..., Ag — 2 (because the minimum number of edges of the genus 
g shape with one face is 2g, and the minimum nonzero P is therefore 4 g). We thus have that, for even 

P, 

n 9 ap)= (e^p*) n ( p2 - ^) 2 )’ 

V t=0 / k= 1 

and it vanishes for odd P, so, after a lower-triangular change of variables {qf ^} —> {b\ 9 ^} and intro¬ 
duction of a normalization factor, its Laplace transform in formula ()3.26[) takes the form 
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(3.27) 


We then use that 


d 1 2 
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or, in the general form, 


( d 2 


- 


\d\ 2 


— (m)‘ 
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m > 1. 


(3.28) 


Now consecutively acting by quadratic differentials in the product, we come to the general repre¬ 
sentation for the one-loop mean, 
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In Sec. [2 we identify the coefficients b^ with ancestor invariants of a cohomological field theory. 
Remark 3. The quantities Cg tS (z) [6] are defined to be the means from (12.611 weighted by z^ ki ^ 2 , 
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From this form and the resolvent representation (|2.13l) it follows immediately that C g \ (z) is related 

(3.29) 


to the one-loop mean wj 9 \e x + e A ) by the simple formula 


w[ 9 \e x + e~ x ) = (e A + e - A )- 1 C 9 ,i((e A + e" A )- 2 ). 

The statement of Corollary 1.5 of |6j is that 

9-1 

C g ,i(z) = P g (z)( 1 - 4 z )-* 9+l/2 where P g (z) = z 2g P^z k , 

k =0 

and Pj. 9 ' 1 are positive integers. The explicit relation between pjf^ and b^ is as follows. If we substitute 
z = (e A + e -A )~ 2 into the above relation and use (13.291) . we obtain that 


Wi 9 \e x + e~ x ) = E P # 
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(3.30) 


A more general relation, which follows from Lemma [6] upon equating all variables A j, reads 
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(3.31) 


and generating functions Cg tS (z) for s-backbone ( 7 -genus chord diagrams are thus determined by g+s —1 
coefficients Pjf’ s ' > - 

Turning back to the one-backbone case, note that, because (e A + e “ A ) 2 = (e A — e “ A ) 2 + 4, we can 
express Pjf^ through and vice versa: 


9-1 

l( 9 ) _ \ ' p( 9 ) /tl—k 

°g-l-k ~ l 

l=k 

A”=E^i-A- 4 )' 

l=k 


ir 


—k 


1 

k ) ’ 


and, in particular, the integrality of the Pjf' 1 proved in [Bj implies that of and vice versa. We prove 
the positivity and integrality of l) 9 ' > in Lemma [8] below; this does not however imply the positivity of 
Pjf^ observed and conjectured in [Bj. 














3.4 The Kontsevich—Penner matrix model 


The matrix model (12.121) manifests many remarkable properties: 

1. It is equivalent USED] to the Hermitian matrix model with the potential determined by the 
Miwa change of the variables tf. = |tr (e A + e~ A ) _fc + \5k, 2 ; 

2. In the special times r = 0,1,..., (13.211) . which is related to the discretization of moduli 
spaces, is equal to the product of two Kontsevich matrix models |9], intertwined by a canonical 
transformation of the variables (see Lemma [3]); 

3. It is the generating function for the discrete volumes N gyS (P \,... , P s ); 

4. It is related to the multiloop correlation functions of the Gaussian model; 

5. It is the generating function for the number of clean Belyi functions, or for the corresponding 
Grothendieck dessins d’enfant [3] (see also ED- 


As concerning the second of the above properties, we have the following exact relation. 


Lemma 3. (J9f) The partition function Z [A] 12.121) expressed in the times Tj^( A) 13.21 1) depends 

ii 


only on the even times T^r(A) and satisfies the following exact relation: 


Z[A] = e 


= ^kp[{T±>] = g C(a7V) p -iV-M p J-K[{T+}]+J-K[{T-}] 


(3.32) 


where Pk [{T^,}] is a free energy of the Kontsevich matrix model 13.171) with being the times of 
the KdV hierarchies and A is the operator of a canonical transformation, 
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(3.33) 


Here C(aN ) is a function depending only on aN that ensures that Pkp[{T^}] = 0 for T± n = 0 and 
B 2 k are the Bernoulli numbers generated by t/(e t — 1) = Ylm=o Bmt m /(m!) and T± k are given by 

dm- 


From this canonical transformation we immediately obtain the (ordinary) graph representation for 
the term P 5)S [{T^}] of the expansion of 


-FkP [{?£}] =Y^ N ' 2 ~ 29 a 2 ~ 29 ~ S J r g A{T?n}}- 

9,s 


Lemma 4. The term ■T,j )S [{T.^ l }] of the genus expansion of the KPMM 12.121) is given by a sum 
of the finite set of graphs G g ^ s described below, with each graph contributing the factor also described 
below decided by the order of the automorphism group of the graph. 


• each vertex Vi, i = 1,... ,q, of a graph G 9tS carries the marking ”+” or genus Qi > 0, 
and has Si endpoints of edges incident to it (2g, — 2 + s* > 0, i.e., all vertices are stable); each 
endpoint of an edge carries a nonnegative integer k^, r = 1,..., Sj, where the superscript + or 
— is determined by the marking of the vertex; 























edges can be external legs (ordinary leaves) with kf } > 0 (we let cq > 0 denote the number 
of such legs incident to the ith vertex), half-edges (dilaton leaves) with kf { > 2 (we let bi > 0 
denote the number of such legs incident to the ith vertex), or internal edges incident either 
to two different vertices or to the same vertex (their two endpoints carry in general different 
numbers kf n and kf 2 i ) (we let U denote the number of internal edge endpoints incident to the 
ith vertex); 

each vertex contributes the Kontsevich intersection index 


T kf."' T k ± . 

l.* =;,*/ gi 


with Y)r=i k r,i = - 3 + Si; 

every internal edge with endpoint markings (kf, kf) or (kf , kf) (two endpoints of such an edge 
can be incident to the same vertex) contribute the factor 


B, 


2 (fcj +kf +1) 


2(fcf + kf + 1) (2 kf + l)!(2Af + 1)! 

and every internal edge with endpoint markings (kf, kf) (two endpoints of such an edge can be 
incident only to distinct vertices with different markings + and —) contributes the factor 

- g 2(fc++fc 2 ~+l) 2 2(k t+ k 2+ 1 ) - l 

2 (kf + k-2 + 1 ) (2 kf + 1 )!( 2 / c 2 + 1 )! 

every half-edge with the marking r* > 2 carries the factor — ^ 9r ± +1 ^ ; 


• every external leg with the marking kf i contribute the corresponding time T) b ± ; 

^ r,i 

• Y^i=i(9i + ii/2 — 1) + 1 = g (the total genus g is equal to the sum of internal genera plus the 
number of loops in the graph); 

• l a i = s (th e total number of external legs is fixed and equal to s); 

From the above formulas, we have that 


s \l\ \B\ 

£g* = 3 s -3 + s-£(l+ tj"f + +fc“) - Y.( k f“' ~ !)■ 

j=1 j=l j=1 


(3.34) 


where, disregarding the vertex labels, kf xt > 0 are indices of the external edges, kff > 0 and kff if 0 
are indices of endpoints of the internal edges, > 2 are indices of half-edges, and \L\ and \B\ are 
the cardinalities of the respective sets of internal edges and half-edges of the graph. 

The proof is just another application of Wick’s theorem, now in the form of exponential of a 
linear-quadratic differential operator: when acting on the exponential of any combination of variables, 
the quadratic part of (any) linear-quadratic differential operator defines the pairing; the sum over all 
possible pairings of an exponentiated quantity (which is in our case the sum of free energies of two 
Kontsevich models) is in turn the exponential of the sum over all connected pairings of these quantities. 
The linear part of the differential operator just produces constant shifts of the higher times, which can 
also be described by insertions of half-edges. We present a typical term in the resulting free-energy 
expansion in Fig. |4j 

This lemma immediately implies the corollary 






Figure 4: The typical diagram from the graph expansion G g , s : it contains three vertices with labels (<?i,+), 
{ 92 ,—), (0, +) and two loops, so the total genus is g = g\ + <72 + 0 + 2. The number of external legs is 
two, so the highest possible total degree of T ± is 6g — 6 + 4 = 651 + 632 + 10. The actual total degree is 
651 + 6(72 — 6 and it corresponds to a reduction of level eight. Recall that every dilaton leaf with marking 
k reduces the total degree in times by 2 [fc — 1 ], every insertion of an edge with markings (+,£+ reduces 
the total degree by 2[fci + fc 2 + 1]. We have four edges and one dilaton leaf, the total reduction of degree is 
2[3 - 1] + 2[0 + 0 + 1] + 2[0 + 0 + 1] + 2[0 + 1 + 1] + 2[0 + 1 + 1] = 16, as expected. 


Corollary 2. The quantities P" 9iS [{T^}] are polynomials in the times such that, for every mono- 
mial T 2 ni ■ ■ ■ T 2n s we have that 1 n * < 3g — 3 + s, and the highest term with Y2i=i ni = 3g — 3 + s 
is 

S S 

(r ni ■ ■ ■ T ns ) g (jj T+. + T 2n .) . 

Z—1 i=l 

This also implies that all discrete volumes Ng )S (Pi ,..., P s ) are Z 2 - quasi-polynomials in Pf. 

Proof. The discrete volumes N g ^(P \,..., P s ) depend only on even powers of P/ because P g>s depend 
only on even times T+ the quasi-polynomiality follows immediately from the fact that P+ are 
polynomials in T 2 + 


Remark 4. Note that the quadratic part of the differential operator (|3.33l) matches the alternating 
structure of (|3.23l) . Indeed, the Bernoulli numbers P 2n manifest an alternating-sign structure, P 2n 
are positive for odd n and negative for even n, 


-E>2n — (—1) 
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2(2n)! 
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Every insertion of an edge with a contributing P 2n decreases the total power k\ + • • • + k s of the 
corresponding monomial T .exactly by n. With the overall minus sign in front of A in (13.3311 
taken into account, this gives the desired alternating-sign structure of (13.2311 up to constant shifts of 
the higher times 2+ with k > 2 ensured by the dilaton leaves. 


3.5 The times for the multi-resolvents 

We now consider the general s-resolvent case. We reformulate Lemma [4] as the following statement 














Lemma 5. From (12.141) we have that the (stable) loop means (with 2g + s — 2 > 1) are polynomials 
in the times for any genus g and number of loops s: 


W^(e Al + e“ A \... ,e A ‘ + = F g , s ({tf .(Xj)}), 


where we have to resubstitute the times T^ d for the new times 
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e Ai ± 1 


(3.35) 


(3.36) 


which are derivatives of the times T^. Note that all the times t± d+1 (X) are now strictly skew-symmetric 
with respect to the change of variables A —>• — A. 


The two nonstable loop means are 
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The proof is the observation that the stable loop means are related to the terms of the free-energy 
expansion of the KPMM by (12.141) and the action of the differential operators in this formula makes 
the time change (|3.36l) . The first nonstable loop mean comes just from the Catalan number 

) is the 


counting and the second one, w!f\ comes from the fact that ^ 0 , 2 (Ai, A 2 ) = — log(l — e 


—Ai—A2 


term from the normalisation factor in the KPMM and acting by Ai - Ai Ai) — on it we 
obtain (13.381) . 

We obtain another useful set of times exploiting formula (13.271) from Example [T] First, we obviusly 
have that 
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with some integer coefficients q 3 .d and q 3 .d- The relation (13.401) follows from 
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g A — 1 e A + 1 e A — e A 
and from another useful representation: 
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We can therefore equivalently expand F gyS ({tf n . +1 (Xj)}) in the variables 


(c^ + 6 

s k,fiW ■— ^ g A _ g— A^2fc+3 ’ ^ — 0, . . . , 3g + S — 3, /A 0, 1. 
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In the next section we demonstrate that the coefficients of these expansions are related to the ancestor 
invariants of a CohFT. 


We now present the general structure of the stable multiloop mean. 































Lemma 6. The general expression for a stable (2g+s—3 > 0) loop mean wj^ (e Al +e Al ,..., e As + 
e _As ) in terms of the variables A) given by (13.421) reads: 

S 

W^(e^ + e _Al ,..., e A ‘ + e~ x -) = Sk^j), (3-43) 

k,0 ' 3=1 

where kj and j3j are subject to the restrictions: 

J s s s 

2 9 ~ 1 + 2 kj < 3g + s — 3, Pj = 0 mod 2. (3.44) 

3 = 1 4=1 4=1 

The proof for the restrictions follows from two considerations: first, as A j —>• oo uniformly for all 
j, that is, we scale A j —>• A j + R, every edge contributes a factor e~ 2R plus s factors e~ R due to 
the derivatives. The minimum number of edges (for a shape with one vertex) is 2g + s — 1, so the 
minimum factor appearing is e (- 4 9~ 3s + 2 ) R whereas Sfc i( g(A) scale as e (^ 3 ~ 2k+ ^ R , which results in the 
lower estimate. The upper estimate follows from the analysis of the pole behaviour as A j = 0. On the 
one hand, Sk t p(\) ~ ^-2^-3 as 0 irrespectively of /?; on the other hand, from the relation to the 
Kontsevich model we can conclude that the pole structure of the derivatives of the Kontsevich KdV 
times is t c j :j (Xj) ~ X rj 3 with dj < 3g + s — 3 and therefore )TT dj = kj, which leads to the 
upper estimate. That the sum of the /5j ’s is even follows from the fact that the total expression is 
symmetric with respect to the total change of the times T± —» T R ; under this change, the variables 
Sfc i( g(A) behave as Sk,p( A) —>• (—l) /3 s/ Cji a(A), so the sum of the beta factors must be even. 


4 Cohomological field theory from discrete volumes 

In this section we describe a cohomological field theory (CohFT) associated to the discrete volumes. 
A dimension d Frobenius manifold structure is equivalent to a CohFT for a dimension d vector space 
H with basis {e a } and metric r/. The main result is Theorem [2] which shows that the quasi-polynomial 
discrete volumes are equivalent to the correlators in the CohFT associated to the Hurwitz Frobenius 
manifold iL 0 ,(i,i) described in the introduction. We give two proofs of the genus 0 case of this, because 
one of the proofs is constructive and the other generalises to all genus. The constructive proof gives an 
explicit realisation of the genus 0 CohFT, and it proves that the CohFT satisfies good properties—in 
particular it is a homogenous CohFT. The homogeneity condition goes some way towards explaining 
the appearance of x{-Mg,n) as the primary correlators in the CohFT. 

4.1 Cohomological field theories 

Given a complex vector space H equipped with a complex metric 77 , a CohFT is a sequence of S s - 
equivariant linear maps 

4.1^' -+H*(Mg, a ) 

which satisfy the following compatibility conditions with respect to inclusion of strata. Any partition 
into two disjoint subsets I U J = {1,... , s} defines a map 

■ -M gu \I\ + l X A4 92j |J| + 1 ->• Mg, s 


and 


| (^) Vi (X) Vj 

\iei j&J 


4>llg,s( v 1 ® ® v s) = J 9i,|/|+1 ® 4,,|J|+l 



where A = Y^ a a p a ^e a < 8 > with respect to a basis {e a } of H. The map 

V’ . J^Ag— l,s+2 ^ ^Ag,s 

induces 

^*Ig,s ( v i ® • • • <8> v s ) = 4_i jS+2 (vi <S> ■ ■ ■ <S> v s <S> A). 

The three-point function /0 3 together with the metric 77 induces a product • on H: 

u»v = ^2 -^ 0,3 e a )rf a ^ep 

a,(3 

where /o^ takes its values in C. A vector eo satisfying 

to, 3 (vi < 8 > u 2 < 8 ) e 0 ) = 77 ( 1*1 ( 8 ) u 2 ), Vui, u 2 € H 


is the identity element for the product on H. 

One further condition, which is not an axiom for a general CohFT, is satisfied by the CohFT 
we consider here (and more generally, for example it is satisfied by Gromov-Witten invariants). For 
s > 3, the forgetful map 

7T : Mg t s+l A4g,s 

induces 

Ig,s+ i(vi < 8 ) • • • < 8 > v s ( 8 ) e 0 ) = TT*Ig )S {yi < 8 > ■ ■ • < 8 > v s ). 


4.2 Quasipolynomials and ancestor invariants 

The discrete volumes N gjS (Pi, ..., P s ) are mod 2 even quasi-polynomials, i.e. it is an even polynomial 
on each coset of 2Z S C Z s . Define a basis of mod 2 even quasi-polynomials induced (via tensor 
product) from the following single-variable basis Pk,a{b) for k = 0,1, 2,... and a = 0,1. 


Po,o(b) 


1, b even 
0 , b odd 


Po,i(b) 


0 , b even 
1, b odd 




The theorem can be proven as an application of m which identifies theories with a special class 
of spectral curves with semisimple CohFTs. We need to check that the spectral curve satisfies the 
necessary conditions. The outcome of applying m is non-constructive so we prove the genus 0 case 
in a different way that enables an explicit realisation of the CohFT. 

4.3 Genus 0 reconstruction. 

The primary correlators of a CohFT are defined by: 

Y,., ■= [_ 4« : H®* C 

J -Mg,s 

which assemble into the generating function 

F(t 0 ,...,tr>-i) = J2 n2 ~ 29 ^ Y 3,s = Y j n2 ~ 29F 9 

where (to, •••, tn-i) hi H* is the dual basis of {eo, eo-i}- The genus 0 part To is the prepotential 
of the CohFT. 

The genus zero primary correlators satisfy relations coming from Keel’s relations in H*(Ai q jS ) [29] 
which we will not write here, and instead use the fact that they are equivalent to the prepotential Fq 
satisfying the WDVV equation. When dim H < 3 (as is the case here) the WDVV equation is trivial. 
Conversely, Manin defines a system of abstract correlation functions to be a collection of symmetric 
functions Yq s : H® s —>• C that form a prepotential Fq satisfying the WDVV equation. 

Theorem 3 Manin [33] Theorem III. 4.3. One can uniquely reconstruct a genus 0 CohFT from 
abstract correlation functions. 

The Deligne-Mumford compactification M g , s possesses a natural stratification indexed by dual graphs. 
The dual graph of S £ Aig )S has vertices corresponding to the irreducible components of E and assigned 
genus, edges corresponding to the nodes of E, and a tail —an edge with an open end (no vertex)— 
corresponding to each labeled point of E. If T is a dual graph of type (g,s), then the collection of 
curves Dp whose associated dual graph is T forms a stratum of M g , s - The closure D p = Up/ < pDp/, 
where the partial ordering is given by edge contraction, represents an element of H*(Ai g s ). Keel [29] 
proved that H*(M. 0 s ) is generated by D p and gave all relations. 

The proof of Theorem [3] uses that 

[_ Io,s( v i ®---®v s )= (5?) Yq\ v \ ( 6?) Vi <g> A s|Er| 

vGV r \i= 1 

which defines evaluation of a cohomology class on boundary strata tautologically from the definition 
of a CohFT. Since H*(M.q iS ) is generated by its boundary strata, and relations in H*(M. 0 s ) agree 
with the relations satisfied by abstract correlation functions, this is enough to prove the theorem. 

In particular, the primary invariants 

To, 3 (eo e 0 (8) ei) = 1 = V 0 , 3 (ei <8 ei (8) ei), V 0)S (e 0 0 anything) = 0, s>3 

Yo, s (ef s ) = N 0>s (0,...,0) = x(M 0>s ) s > 3 (4.46) 

define a genus 0 CohFT. The claim of Theorem [2] that the coefficients of Nq )S are the correlators of a 
CohFT, is an overdetermined system, since the CohFT has been uniquely determined by the constant 
terms Nq jS (0, ..., 0). To prove that the coefficients of Nq are indeed the correlators of this CohFT we 
use homogeneity of the correlation functions and the CohFT. 



4.4 A homogeneous CohFT 


A CohFT is conformal if its prepotential is quasihomogeneous with respect to a vector field E known 
as the Euler vector field: 

E ■ F 0 = (3 - d)F 0 + Q(t) (4.47) 

where Q is a quadratic polynomial in t = (to, to-i)- Using the genus 0 reconstruction in Theorem[3l 
Manin proved that a conformal CohFT induces the following push-forward condition on the genus 0 
CohFT. 

Let £ be any vector field on H which we consider to be a manifold with coordinates to, to—i € H*. 
The Lie derivative with respect to £ of the correlators of a CohFT I gtS induces a natural action on the 
CohFT: 

S 

■ -®V S ) = deg s (ui <g) • • •®v s )-'y' j Ig,s(vi®»®[€,Vj\<8>»®Vs)+n*Ig,s+i(vi® — ®Vs®€) 

3 = 1 


where 7r : M gtS+ i —>• M.g )S is the forgetful map. 

Here I is a (H*(Ai o i < j )-valued) tensor on H which is acted on infinitesimally by the vector field £ 
on H. 

A CohFT is homogeneous of weight d if 

(E.I) g , s = ((g-l)d + s)I aiS (4.48) 

If the preprotential satisfies the homogeneity condition (|4.47|) then the proof of Theorem [3] can be 
used to prove that the genus 0 CohFT is homogeneous. The Lie derivative of the bivector A dual to 
the metric r] on H can be calculated in flat coordinates 

Ce ■ A = Ce ■ rf^ei ® ej = rf j ([E , e*] <g> ej + e* <g> [E, efi) = (d — 2)r/*- 5 e i ® ej = (d — 2) A 

where we have used a choice of flat coordinates 1141 with respect to which rj = Si D—i—i an< 4 

E = ^2{aiU + 


where otj + otD-i —i = 2 — d. The common weight d in the homogeneous conditions (14.471) and (14.48[) 
arises from the following compatibility for pull-backs in a homogeneous CoHFT. 


p*E ■ I, 


9,s 


E ' P Ig,s = E ■ {Ig lt si- 1-1 ® ■^J2,S2+l(^)) 

{E • -^)g 1 ,s 1 +l ® -^92,32+1 (^) 4" -^91,51+1 ® {E ■ L)g 2 ,S 2 + 1 (A) + ® I g2tS2 + l(£,E ' A) 

[((^l — 1)^ + Si + 1 ) + ((<72 — 1 )d + S2 + 1 ) + (d — 2)\Ig ljSl+ i ® Ig 2 ,s 2 + l(A) 

(id — E)d + s)I gitSl + 1 <8> Ig 2 ,s 2 - 1-1 (A) 


4.5 Proof of Theorem [2] in genus 0. 

We are now in a position to prove the genus 0 case of Theorem [21 The idea is to produce a prepotential 
from the primary terms of iVo iS (Pi,... , P s ), i.e., from their constant terms, which uniquely (and 
constructively) determines a genus 0 CohFT. Moreover, quasihomogeneity of the prepotential produces 
a homogeneous CohFT. The higher coefficients of Nq <s (Pi, ..., P s ) satisfy a homogeneity condition that 
forces them to be the correlation functions of the homogeneous CohFT. 








The prepotential 


To — ^ -yko.s — 9 ^ 0 ^! + ^yT/o,s(0)if — -tlh + -(1 + ti) 2 log(l + ii) — ~h — -t 2 (4.49) 


s>3 
\s—3/ 


assembled from IVo iS (0) = (—l) s J (s — 3)! is quasihomogeneous with respect to the Euler vector held 
E = t omE + 2 ( 1 + t i)mi : 

E ■ Fq = AFq + t\ + fg. 

This ensures that the genus 0 CohFT /o jS produced from Theorem [3] satisfies 

K*I g ,s+i(es (8) ei) = i(l - g + s - deg - ^ a ik )Ig }S (e s ) (4.50) 

where es = eq <8>... (8) e^. and ao = 1, a± = 2 are the coefficients of E. 

The CohFT also satishes the pull-back condition described above 


I g ,s+ i(ffi ® ••• ® v s (8) e 0 ) = vr*/ S!S (ui <8> ... <8> v s ). 

It is proved in genus 0 via construction, using the fact that the only non-zero genus 0 primary invariant 
with eo as an input is /o, 3 (eo <8> eo <8> ei) = 1. 

Teleman [33] produces a unique homogeneous CohFT for all genus extending the homogeneous 
genus 0 theory. 


Theorem 4 Teleman [44]. A semi-simple homogenous CohFT with flat identity is uniquely and 
explicitly reconstructible from genus zero data. 


Thus, given the genus 0 primary invariants jVo,. s (0) there is a unique homogenous CohFT with hat 
identity. Below we will see that its correlators agree with the coefficients of N g s (P \,... ,P S ). 

The pushforward relation (14.50ft expressed in terms of correlators is 


IM 


s „ s 

Ig,s+i(es ® ei) = _ 7T*( Y[iF*^i) ki P^Dj ■ ]^7rV^“Mlg )H _i(es <8> e{) 

7,s+l i —1 J M.g,s l j 

p s s 

= _ \\^^Jg,s+i{es® ei) + • Ig iS+ i)(e s <8> ei) 

J M g ,s i 


+ Xg,. 


i= 1 


Mn 


J g,s(es) n — 4* ( e s\{j} ® e *i ) n ^ 


1=1 


j =1 J Mg,‘ 


i= 1 


which uses 

1>i = + Di=> $ = (7 r*^i) k P A ■ £ Vi) fc_1_m 

and 

Dj ■ Ig, 3 +i(e s <8> ei) = I g<s <8> I 0 ,3(e S \{j } <8> A <8> e*. <8> ei) = T s , s (e 5 \{ j } <8> e*) 

for e* = ei_j. 

The condition E- Fq = 4Eo + + fg on A r o,, s (0) is a specialisation to g = 0 and Ej = 0 of the divisor 

equation [39] 


N g , s+1 (0,P 1 ,...,P s ) = J2J2 k Ng, s (P 1 ,...,P s )\p j=k + ( l^Pj + Xg,,) Ng, s (P 1 ,...,P s ) (4.51) 


j =i fc=i 


l=i 




which is exactly the pushforward relation on correlators (14.511) . 

The flat identity pull-back condition is known as the string equation on correlators for 2g— 2+s > 0: 


■' -Mg } 3 + 1 

and agrees with the recursion 


I g ,s+i(vi ® ® ® e 0 ) JJ ^ 


2=1 



Ig,s( v 1 ® ® V s ) V’f’ 5t ' J 

2=1 


N, 


g,s+l 


(1 ,Pl, ■■ ■ ,P S ) = kN g iS {Pl, ■ • • ,P Si 


P,=fc 


(4.52) 


j= ifc=l 

In particular, this proves the genus 0 case of Theorem [2] since the recursions (|4.5ip and (14.521) uniquely 
determine the correlators of Iq iS and A r o.. s (Pi ,..., P s ). 


4.6 Explicit description of CohFT in genus 0 


This constructive proof gives rise to a rather explicit description of the genus 0 classes /o,s(es) £ 
H*(M. o iS ). It is sufficient to describe the pairing of the class h). s { e s) with all strata Dr since this 
uniquely characterises /o,s(es). 


Proposition 1. We have that 



(4.53) 


where evaluation of x(Z?r) or 0 can be characterised explicitly in terms of T and e$ = eq <g> ... <8) e, s . 


Proof. If T consists of a single vertex, so that Dr = Mo,s then (|4.53p holds by construction—it is the 
constant term of Ntr s (P \...., P s ). In other words 



x(Mo, s ) e s = ef s 
0 otherwise 


Given a dual graph T of type (0, s), its associated strata is a product of (open) moduli spaces, one 
factor for each vertex of T, 

D r = A4o, Sl x ... x A4o, S | V(r) | 
where Sj is the valence of the jth vertex of T. Hence 

x(Dr) = X(M 0 , S1 ) ■ • x(Mo, s J. 

In terms of the primary correlators Yo^ Si = J-j^ Q Io )Si ■ H® Si —> C 

[_ IoAes) = ® Y oM (e s ® = c(T, e s ) • x(M 0 , sl ) ■ ... ■ x(M 0 , Sk ) = c(T, e s ) ■ x(D r ) 

jDr v&V v 


for some c(T, e^) E Z since each primary correlator is either zero or an Euler characteristic. Given es, 
if there is an assignment of eo <S> ei or ei <S> eo to each interior edge which gives a non-zero evaluation, 
then it is unique since T is a tree. Hence c(T, es) = 1 or 0 as required. □ 


Remark 5. More generally the same argument shows that for any genus 

[_ IgA e s ) = C (T, es)x(-Dr) 

J D r 

for c(T, eg) G Z-°. The integer c(T,e 5 ) can take values other than 0 or 1 since there can be more 
than one assignment of eg (g> ei or ei <8> eo to each interior edge which gives a non-zero evaluation. 










4.7 General proof of Theorem [2] using DOSS method |17| . 


In higher genus we have agreement of relations among the correlators of the CohFT evaluated on 
boundary classes and relations among the coefficients of N gs (P ],..., P s ), but this is not quite enough 
to prove that they coincide. Instead we apply the results of m where it is shown that if the spec¬ 
tral curve satisfies a constraint then the Givental reconstruction of higher genus correlators can be 
understood in terms of graphs, and the same graphs are used to calculate topological recursion. 


Dunin-Barkowsky, Orantin, Shadrin and Spitz m using Eynard’s technique of [T9] associated to 
any semi-simple CohFT a local spectral curve (£, B, x, y ). The i?-matrix gives rise to the bidifferential 
B on the spectral curve 


B % ' 3 = 

P,q 


\z p w q 1 




z + w 


(4.54) 


where Bp? q are coefficients of an asymptotic expansion of the Laplace transform of the regular part of B 
expressed in terms of the local coordinates Sj = \Jx — x{ai) where dx(a{) = 0. See [IT] for details. The 
/2-matrix together with the transition matrix T from a flat to a normalised canonical bases gives rise 
to the meromorphic differential ydx in terms of s j. In particular, this implies a compatibility condition 
between the differential ydx and the differential B. It is given by (|4.55|) below. Most spectral curves 
will fail this condition. 


One can apply m in either direction, beginning with a semi-simple CohFT or a spectral curve. 
The prepotential Fq (|4.49p gives rise to a semi-simple CohFT. It can be used to produce the //-matrix 
and transition matrix \k and hence the spectral curve. Instead, since we already have a candidate 
for the spectral curve we will start with the spectral curve and apply m to get the coefficients of 
N g s {Pu • • •, Ps ) as ancestor invariants of a CohFT. Since it agrees with the CohFT above in genus 0, 
by uniqueness it is the same CohFT produced by Teleman’s theorem. 

The spectral curves for the discrete volumes and Gromov-Witten invariants of P 1 are similar: 


x = z + 1/z , y = z, 


B = 


dzdz' 
(z — z') 2 


x = z + l/z, y = \ogz, 


B = 


dzdz' 
(z — z') 2 


and since x and B determine the //-matrix, it is the same for both curves. The //-matrix for Gromov- 
Witten invariants of P 1 was calculated explicitly in [IT]: 


R(u) = ^R k u k , R k 

k =0 


(2k - 1)!!(2 k - 3)!! 
¥^k\ 


-1 (—l) fc+1 2 ki \ 

2 ki (-l) fc+1 J ' 


The results of Em can be applied to spectral curves such that a Laplace transform of ydx is related 
to this R -matrix (which is essentially the Laplace transform of the regular part of the bidifferential). 

For local coordinates s*, i = 1,2 near x = ±2 given by x = sf ± 2 


1 OO 

v = l + s i + t; s i + y~1(—l) fc 1 


k =1 


(2k — 3)!! 2fc+1 
2 3k k\ bl 


V 


— 1 + is 2 + -s 


2 

2 


OO 




(2k — 3)!! n 2k+i 

2 3k k\ 2 


v T' 1 k =0 















(fc+1) 


{ydx) 2 = 


2^ 


/ OO 
e ~ u (x + 2)ydx —i ^ . 

' 2 k = 0 


(2fc + l)!!(2ife-3)!! 
2 ik+l k\ U 


We use the convention (—1)!! = 1, (—3)!! = —1 and ~ means Poincare asymptotic in the parameter u. 

The compatibility condition between the differential ydx and the differential (appearing here in 
terms of the i?-matrix which is related to the bidifferential via the Laplace transform) is 


1 

71 


i 


i 


■ Tam 


which uses the first row of the transition matrix 


(ydx) i (ydx)r 



1 

1 


i 

—i 


(4.55) 


One can check directly that it is satisfied for x = z + 1/z, y = z, B = dzdz'/(z — z') 2 . 
From this, m supplies the times 


and the main result 



Wj 3) (x i,...,s a ) = XMsII&i. 0 * 


k,a 


2—1 


(4.56) 


where the coefficients are the ancestor invariants (14.451) . As remarked above, the CohFT produced 
this way necessarily coincides with the homogeneous CohFT produced by Teleman’s theorem since 
they both use Givental reconstruction and the same initial data. 


4.8 Euler characteristic of the moduli space of curves 


The primary correlators of the CohFT described above are given by the orbifold Euler characteristic 
x(-M. g ,s)- This is perhaps not so deep in genus 0 where it is simply the statement that 



d s F 

~dt[ 


= x(Mo, s ) 

(to,ti)=(0,l) 


for F = \tl oto,i + 7o,i log to,l given in (11.11) . It is a deeper property that 


[_ I g , s (eT) = x(M g , s ) 

and it is the purpose of this section to show that homogeneity of the CohFT goes some way towards 
explaining the appearance of x(A 4 3iS ). 

The exact sequence of mapping class groups 












implies x(r®) = x(^ S g + )/x(C ~ {Pu •■-,Ps})- Since the (orbifold) Euler characteristic is x(^g.s) = 
X(rp then 

x(Mg iS+ i) = (2 - 2g - s) ■ x(M g>s ) for 2g - 2 + s > 0 (4.57) 

and x(M gtl ) = C(1 - 2g) mM- 

The relation (14.57|) arises naturally out of the push-forward relation satisfied by the CohFT which 
restricts to the top degree terms. 

Tr*Ig, s+ i(e h ® ... ® e is ® ei) = i (l -g + s- deg - ^ a ik ) I g A e h <® ■ ■ • ® e is ) (4.58) 

Note that it is usually more natural to consider the degree 0 part of the CohFT —a topological field 
theory— since the pull-back relations satisfied by a CohFT restrict to degree 0. 

Denote the degree 3g — 3 + s part of I g>s (e f s ) by c 9)S so 

Ig,s{ e V) = e (d + s) 2 9 + ... + c g , s {pt} 

where e(g + s) € {0,1} is the mod 2 reduction of g + s. Note that all insertions are necessarily e\ since 
any eo insertion comes from a pull-back 

deg/ 5jS+ i(e 0 ® es) = degn* I g . s (e s ) = deg7r* I g>s (es) <3sr-3 + s + l 

and hence is not of top degree. The top degree class is well-behaved under push-forward. 

tt*c s , s +i{pt} = ^ (1 - g + s - (3g - 3 + s) - 2s) c g , s {pt} = (2 - 2 g - s)c 5iS {pt} 


hence 

Cg,s+i — (2 2 g s^Cg tS =r- c g s — Cg ■ x(A4 g,s)- 

This explains the appearance of the Euler characteristic x(A4 ffiS ), up to a factor, in this CohFT. It 
turns out that c g = 1 but this needs the discrete volumes to prove. 


4.9 Ancestor invariants and Gaussian means 

From Lemma [2] and the formulae (|3.27l) and (|3.41l) . we immediately obtain formula for the loop means 
in terms of the ancestor invariants. 


Theorem 5. We have the following explicit relation between the ancestor invariants 114.451) of a 
CohFT and the Gaussian means: 


where 


wA\e Xl + e" Al ,...,e As + e" As ) = E c |, fi II 

k,a 

2 1_2r (2r + l)s ri0 (A), k = 2 r, a = 0; 

2 _2r (2r + l)s r ,i(A), k = 2r, a = 1; 

A) — 

2 -^(2r + 2)(2r + 3)s r+ i,i(A), k = 2r + l, a = 0; 

2 _2r_1 s r! o(A), fc = 2r + l, a = 1, 


(4.59) 


(4.60) 


and s Tt p( A), /3 = 0,1, are defined in 113.421) . 













Example 2. The topological (= degree zero) part of the CohFT is 


Ig,s(ea i <8> ■ ■ ■ <8> e Qs ) = e(a) 2 9 + higher degree terms 

S 

where e(a) = a* (mod 2) is 0 or 1. This explains the asymptotic behaviour of the topological 


2=1 


invariants W^ at their poles. 

Example 3. If {eo,ei} is a basis of H corresponding to flat coordinates then 

[_ 4,(ef) = X (M g , s ). 

J A4.g } s 

This uses the fact that N gtS (0, 0,..., 0) = x{Mg,s) and 

' 1, (k,a) = (0,0) 


Pk,a( 0 ) = 


0, otherwise. 


We thus identify the coefficients 6| - of the expansions (13.431) with (linear combinations) of the 

ancestor invariants cP _ using the identification (14.601) : for s = 1, we have 
k,a 


b% = 2 1 2r (2r + l)c? 2rfi + 2 1 2r cf r+11 , 

b 9 r , i = 2 _2r (2r + 1)4. x + 2 2 ~ 2r 2r(2r + 1)4_ 1>0 . 


5 The one-backbone case. The Harer—Zagier formula 

5.1 Times in the one-backbone case. 

The formulas (|3.35l) and (|3.36l) become especially useful in the one-backbone case. Then, we have just 
one variable A and just one set of times, and we always have the sum (13.391) in the r.h.s. of (13.351) . It 
is useful to introduce the special notation for the one-loop resolvent 

L g (x) := W[ 9 \x). (5.61) 

We have just one external edge in this case, so from ()3.35jl we merely have 

L„(e* + e-*)=f>l Y~p 

r =0 

where x g .\,r are (conjecturally positive) rational numbers, Xggy = {Tz g - 2 ) g - 

5.2 The recurrence relations 

The expansion coefficients c g (k ) (in the original notation) of the loop mean 

OO OO 

L g (w ) = ^ yr a_1 (trH 2k ^ = ^ c g (k)w^ 2k ^ 1 (5.63) 

k=2g 9 k=2 g 

satisfy the celebrated Harer-Zagier three-term recursion relation [27] : 

(n + 1 )c g (n) = 2(2 n — l)c 5 (n — 1) + (2 n — l)(n — l)(2n — 3)c ff _i(n — 2). 


K gh,r 


~(d — r)! e A — e 


-A 


d_ 

d\ 


za-zr -|-i 


s 2A_i’ 


d = 3g-2, 


(5.62) 


(5.64) 


















Lemma 7 . The generating function L g (w) satisfy the differential recurrence relation 

L'g{w) = (w 2 - 4 )L' g+1 (w) - wL g+ i(w), (5.65) 

which, upon substitution w = e A + e _A , takes the form 


L"'{er + e _A ) = (e A - e^L' 1 (e A + e“ A ) - (e A + e~ x )L g+1 (e x + e" A ) 


(5.66) 


Note the appearance of the factor (e A — e A ) 2 of the first term in the r.h.s. It simplifies the further 
recurrence procedure drastically. 

The next proposition is a corollary of Lemma [6J 


Proposition 2. The general form of L g (e A + e A ) is (cf. Example QP 


kg (e A + e A ) — 


9-1 

-E 




(e A — e a ) 4 9+! E-j' (e A — e A ) 2fc 


(5.67) 


Remark 6. Note that the coefficients l)f ' 1 are the coefficients n g {n ) appeared in the original pa¬ 
per m by Harer and Zagier. But the three-term recurrence relation on them was not written there. 
We present it below. This relation was also found in 


Proposition 3. The coefficients b jr 9 * in Proposition^ satisfy the three-term recurrence relation: 
(4 g + 2 k + 6 )b { k 9+1) = (4 g + 2k + 1)(4 5 + 2 k + 3) [(4 g + 2k + 2)b[ 9) + 4(4 g + 2k- l)&jf_M . (5.68) 


Proof. The proof is the direct substitution, we just mention the basic steps. We first find the 
action of the third derivative on any monomial (e A — e -A ) -< : 


d 


e A _ e -A ( J\ 


(e A — e~ x Y 


— — t(t + 2) 


(t + 1) 


e A + e A 
(e A — e _A )*+ 4 


+ 4 (t + 4) 


e A + e A 
(e A — e _A ) t+6 


The result thus contains two terms, whereas the action of the operator in the right-hand side of (15.6611 
on the same monomial is merely 


d 


(e A -e- A )--(e A + e 


-A', 


1 


(e A — e A )* 


— ~ ( i + 1) 


e A + e A 
(e A — e~ x Y 


and it contains just one term. Substituting expansions (15.6711 in both sides of (15.66|) and equating 
coefficients at equal powers of , e \^ e _x\ s we obtain (15.6811 . which completes the proof. 

Lemma 8. All the coefficients b^ are positive integers. 

Proof. The positivity obviously follows from the fact that the both coefficients are positive in the 
r.h.s. of the recurrence relation (|5.68ll . The integrality follows from the comparison with the expansion 
of the original Harer-Zagier function L g (e x + e _A ) at large positive A. Indeed, comparing formulas 
(|5.63H and ([5.6711 and expanding in e~ x it is straightforward to obtain that 


oo oo 


L g (e x + e x ) = 

k=2g p =0 
9“1 oo 

- EE4 S) 

k =0 p =0 


P 

2k + p 


0 -(2k+2p+l)\ 


P 

Ag + 2k + p 


o —( 4 9+2fc-t-2p-|-l)A 


(5.69) 

































All the Harer-Zagier coefficients c g (k) are positive integers. If we use formula (J5.69I) to express b^ 
through c g (k) equating the terms at equal powers of e _A , we obtain for the coefficient of the 

equality 

b o 9) = c q( 2 9) 


and b^ is therefore an integer. The equality at e ( 4 9+ 3 )- >i reads 




1 

4g + 1 


Cg( 2 ff), 


so l)f' 1 is an integer as well, etc. Every first appearance of any coefficient bjr 4 in this chain of recurrence 
equations is with the unit coefficient, which enables us to express b^ as a polynomial in the preceding 

b p 9 \ P = 0, • • • j k — 1, and in c g {2g + q), q = 0,..., k, with integer coefficients, so, obviously, b^ is an 
integer, which completes the proof. 

Conjecture 1. Because we have identified the coefficients b^\ with combinations of the ancestor 

J k,p 

invariants in the general (multi-backbone) case as well, we put forward the conjecture that all these 
coefficients are positive integers. 

In order to simplify expressions it is convenient to introduce the renormalized bjf^ by the formula 

b[ 9) = bjf } (6 g -1-2jfe)!!. (5.70) 

In terms of these renormalized quantities, the general recursion relation ()5.68|) looks especially simple: 


(2 g + k + 3 )b<f +1) = (2 g + k + l)bjf + 2b[ 9 > v b^ 


,(ff) 


,(s) 




(5.71) 


5.3 Examples 


Note first that b^ = 1, which follows from the simple diagrammatic calculation in the genus-one case. 
We obtain the coefficients b^ and 6 ff _ 1 from the simple one-term recurrence relations obtained by 
setting respectively k = 0 and k = g in (|5.68l) : 


(4g + 6)b { 0 9+1) = (4 5 -l)(4 5 + 3)(4 5 + 2)6^, 


(6g + 6)b<f + V = 4(6 5 + l)(6 5 + 3)(6 5 -l)6^ 1 


which immediately give 


(9) 2 9 ~ l (6g — 3)!! 

9-1 39 c/! 


b is) = ( 4 ff) ! 

0 89 5 !(2 5 + l)!!' 


Note that in terms of the renormalized quantities (15.701) we merely have 

2 9 - 1 . 1 


b (9) - 
Vi - 


3 9 g\ 


b (9) = 

0 2g + V 


(5.72) 

(5.73) 

(5.74) 

(5.75) 


Now substitute b\ 9 }_ { into (15.621) and evaluate the term without reduction (r = 0) which contains 
the highest Kontsevich coefficient (r^g-^ . The easiest way to do it is to compare the highest-order 
poles in the both expressions. The result reads 


(735-2)5 - 


1 


'a 2^9 2,9 g\ 

and it is in the complete agreement with the standard formulas. 


(5.76) 

















5.4 Solving recursion equations for b^ J \_ k at fixed k and arbitrary g 
5.4.1 The coefficient b^]_ 2 


The recurrence relation for the renormalized b\ 9 l 2 rea< l s: 


(35 + 2)b|f_ + 1 1) = (3 + 2b% 


bL 1 ] = 0 . 


(5.77) 


We try the anzatz b^ 2 = f(g)b ( g 9 } 1 . Then, since bjf +1) = 2b ( g 9 } 1 /(3g — 3), we obtain for f(g) the 
relation 

(35 + 2 )/(<? + 1)^— = 3 5 + 2 /( 5 ), 

35 + 3 

which can be solved by a polynomial substitution provided we can cancel g +1 in the denominator in the 
right-hand side. So, taking f(g) = 05(5 — 1 ), we obtain merely that (35 + 2)(2a/3) = 3 + 2 a (5 —1). For 
the term linear in 5 we have the identity, whereas the constant term gives a = 9/10. The normalization 
condition at 5 = 1 in (15.771) is satisfied, and we therefore obtain a unique solution to the recursion 
relation (15.77[) : 



y+9 - I)'’ 


(9) 

9-1’ 


i.e., 



1 29~ 2 

5 39-2(5-2)!’ 


or 



1 29 - 2 (65 - 5)!! 
5 3»- 2 (5-2)! 


(5.78) 


that 


Expanding in the times t± d+1 ( A) (see (|5.62|) ). for the first “reduction coefficient” i 1 we obtain 


^s.i.i = 7\ l2 9 2 - 75 + 5]x ff) i i0 , g> 2 . 


(5.79) 


5.4.2 The coefficient b 


(9) 

9-3 


The recurrence relation (|5.68l) for k = g — 2 reads 

(65 + 2)6^_ + 2 1) = (6 g ~ IX 65 - 3) [(65 - 2)bf_ 2 + 4(65 - 5 ) 6 ^]. (5.80) 

Recall that due to (|5.78[) the coefficients b^]_ 2 now themselves satisfy the two-term recursion 

(5.81) 


u ( g + 1) _ 2(65 + 1)(65 - 1)(65 - 3) M 
y —1 - y —2 


3(5-1) 

and we can try to find an anzatz of form (15.781) relating now bg 9 } 3 and b \ 9 }_ 2 . This anzatz must have a 
form of a rational function with polynomial of the second order in the numerator. But it must also 
cancel the factor 65 — 5 in the recurrence relation. Moreover, the polynomial in the numerator must 
be divisible by 5 — 2 to cancel the factor 5 — 1 in the denominator of the recurrence relation (I5.8ip . 
So, we try the anzatz 

(n - 9)(n - t\ 

(5.82) 


1(9) _ (g-2)(g-t) (g) 

y~3 - K 77—T %-2 


65-5 

with unknown k and t. Substituting this anzatz into (15.801) . we obtain the simple polynomial equation 
of second order in 5 : 

^(65 + 2 )n(g - t + 1 ) = (65 - 2 ) + 4^(5 - 2)(5 - t) 
whose highest term always matches. A unique solution is t = 1/2, k = 3/10. But if we substitute this 

/ o\ 

solution into expression (j5.74|) for the term 6 q j we observe the mismatch, which is due to the fact that 

























we are able to add to the obtained solution 6^ 3 of the inhomogeneous equation (15.80|) any solution 
3 of the homogeneous recursion relation 


( 6 g + 2)6^_ + 2 1) = 4(65 - 1)(6 5 - 3)(6 5 - 5)6^, 
so, with the proper normalization for the term bj^ , we obtain b^ 3 as a sum of two terms : 


(5.83) 


= (2g :T 2 r (69 .: 7> " - 72 ! , (eg -here (3e - 2) m ft* 2). 


9~ 3 5 2 39-3(c/-3)! 


10(35-2)!! 


(5.84) 


k =3 


The general structure of solution becomes clear: in order to find 4 we first solve two inhomo¬ 
geneous equations with and in the right-hand sides and then again add the solution of the 
homogeneous equation with the proper coefficient to satisfy the initial condition for 6 q . We now 
formulate the general procedure for finding all coefficients at all 5 and fixed k. 


5.4.3 Finding 


We formulate the general procedure in terms of the renormalized coefficients b^\_ k (|5.70l) . The 
procedure consists of several steps. 


Step 1. Solving the homogeneous equation for enumerated by s = 0,1 


(3 g~s + 3)b^ 1} = 2bJfj s _ l5 or b^ 1} = 


'(9) 


9(5 + 1,«) 9 ' S " 1 ’ 


where 


Then, obviously, 


Q{ 9 ,s) :=-(3g - s), s = 0 , 1 ,.... 


g-s-l 


^g-s- 1 - [ II Q( s + 1 + 1 i' 


i -U 


1=1 


:(«+i) 


(5.85) 

(5.86) 

(5.87) 


Step 2. Solving the chain of inhomogeneous equations for b^ s _ 1 _ fc that are solutions of the 
equations 

5(5 + 1 , s + k)b%+V k = q(g, s + k- l)bf_ s _ k + 
by substituting the anzatz 


(9)9(9 - k + 1, a) 


P(s,k— 1) ( 5 ) 


where 


b (s) _ -g-(fir) 

°g-s-l ~ 0 g-s-l 


and all P( s j-) ( 5 ) are polynomials of order k to be determined (-P( Si 0 )( 5 ) = 1 ). 


It is easy to see, by recursion, that the quantities 


b (9) , 

g — s — k 

P(s,k— 1 ) ( 9 ) 


(5.88) 


(5.89) 

(5.90) 


satisfy the recursion relation 


u(9+l) 


f b (g) 

g+l—s—k 

= q 1 {g + 2-k,s) 

°g-s-k 

P(s,k-i)(g + 1 ) 

P(s,k— 1) ( 5 ) 






















whence, for P( s ,k)(g)> we have the recursion 

(3g + 3-fc-s)P (Sjfc )(g + l) = (3g + 1 - k - s)P( Stk _ 1) (g) + 2q(g - k + l, s)P {Stk) (g), P( S)0) = 1, (5.91) 
or, explicitly, 

(3^ + 3-fe-s)P (Sjfc )(5 + I) = (35 + l-A;-s)P (Sjfe _i)( 5 ) + (35-3A; + 3-s)P( S)fc )(5), P (S)0 ) = 1. (5.92) 

Solving this relation, we can find all the polynomials P( s ,k)(g) starting from the zeroth one. This 
recursion always has a solution because every highest term relation (the coefficient of g k+1 ) is satisfied 
automatically and at each step we have a system of k + 1 inhomogeneous linear equations on k + 1 
coefficients of the polynomial P( s ,k)(9 )■ 

Step 3. In light of (15.871) . we now express b ^ s _ fc _ 1 as 

g—s—l—k 

b^ls-k-i = P(s,k)(g) II q(s + l + l,s)b [ Q +1) . (5.93) 

1=1 


for b 


Step 4. Having all the polynomials P^ s k)(g) for s,k < r determined, we define the general solution 

(9) . 

g — r — l ' 


r g—l—r 

b( g-r-l = ^2 P (s,r-s){g) IJ ?( S + 1 + h s)b^ +1) , 


(5.94) 


s=0 


1=1 


where the coefficients bQ S+1 ^ are to be defined from the triangular system of linear equations originating 
from the fact that b^ = 1/(25 + 1) : 


(1) - 1 w ' +1> - ~ - E F(p,,-p,(»+i)S ,p+1) 


b v ; — - b 
l1 0 — D 0 


2s + 3 


J o 


, s > 0. 


(5.95) 


p=0 


We can therefore derive the coefficients Xg i r for an y 9 for a fixed “degree of reduction” r. 
Remark 7. The first few polynomials P( S! k)(g) are 

P{s,o)(g) = 

P(s,i)(g) = ^[35 - a - 3 ], 

, , 9 o 33 3 1 , 11 39 

P/ s 2 )(q) — —9 - 9 - 9 s H-s H-sH-. 

50 70 25 50 35 175 

So, we can conclude that the quantities P( S) k){9 ) are polynomials of degree k in both 5 and s, and 
their general form is 

P(s,k)(g ) = Y a S } - 9 V ■ ( 5 - 96 ) 

i>0, j>0 

i+5 <fc 

From relation (|5.92p we can read the recurrence relations on the coefficients a[ k j: rewriting (15.92[) in 
the form 


(35 — 3 — k — s) 


E ag’ffe+l+sV- E 

i+j<k i+j<k —1 


(fe—1) i j 
a i,j 9S 3 


- 2 Y ,ps ’ 

i-\-j<k —1 


+ — 2A: ^ ntyjj 


a i,j 9 s J , 


(5.97) 


z+_7</c 


which produces the recurrence relations if we equate coefficients of g l s J : 











• for i + j < fc — 1, 


k-j 

3 E 

p=i 


i — 1 
p 


3 + 0-lfe) £ 

p=i+l 


.(fc) 


P 


(fc) 

“pj - 


fc+l-J 

£ 

p=i+l 


P 


M (fc) o (fc-1) /O 7 \ (fc-1) . (fc-1) 

°Pj-i - 3a i-ij - (3 - k ) a ij + %j -1 


= -2ag- 1) -2fcag ) , 


(5.98) 


where the terms with i — 1 are absent for z = 0 and those with j — 1 are absent for j = 0 ; 


• i + j = k, i > 0 , j > 0 : 



O • (k) /. -.s (fc) 0 (fc—1) . (fc—1) . r) 7 (fc) n 

3m ij ~ u + iK+u-i - 3a i-ij + a i,j-i + 2ka ij = °; 

(5.99) 

0, j = fc: 

(fc) . (fc-l) , 07 (fc) n 

a l,fc-l + °0,fc-l + 2ka o,k ~ 

(5.100) 

II 

O 

ri (fc) 0 (fc—1) (fc) 3^ 

5 fc«fc,o = 3 ®fc— 1 , 0 ’ or a fc,o = 5 fc kl ■ 

(5.101) 


We see that the quantities ag' 1 with z+j = fc satisfy the separate recursion relations (j5.99j) - (|5.101|) : 
we can first solve these relations, determining subsequently all lower degree terms from the relation 

flgjggP - 

5.5 Deriving from the graph representation of Lemma [4| 


We now derive the first subleading coefficient g _ 2 (the first ancestor invariant) of the expansion of 
the one-loop mean L g (e^ + e -A ) f)5.67|i using the decomposition formulas and Lemma 01 

We consider only the part with the times T^ k . The highest term for genus g is (Tzg- 2 )gT^ g _ A 

Following Lemma 01 the first-order correction, or the coefficient of Tg_ 6 , comes only from two 
terms: from the graph with one vertex and one internal edge with endpoint markings ( 0 , 0 ) and from 
the graph with one vertex and one half-edge with marking 2 (see Fig. [5]): the corresponding coefficient 
is then 

B 2 3 

^(r 3 g—3T 0 ro ) 5 —1 + -^-(T 3 g _3T2) 9 , (5.102) 

and we need only to know the corresponding intersection indices. Whereas (r 3g _ 3 roTo) 9 _i = (r 3s _ 5 )^- 1 , 
the intersection index (r 3 g _ 3 T 2 ) s must be calculated separately, and we do it in Lemma IA.9I in Ap¬ 
pendix [575] Using formula (IA. 6 I) and that B 2 = 1/24, we obtain that the coefficient in front of Lg _ 6 
is 

i[12 ff 2 -7 ff + 5] (5.103) 

in full agreement with (15.7911 . 


Conclusion 

Application of topological recursion (TR) to constructing generating functions for cohomological field 
theories is becoming an important issue in contemporary mathematical physics (see, e.g., the recent 
paper [22] where all genus all descendants equivariant Gromov-Witten invariants of P 1 were constructed 
using TR. In this respect, it seems interesting to understand the status of Givental-type decompositions 
in the quantum spectral curve approach. 

















Figure 5: The two diagrams contributing to the first-order correction to the Kontsevich potential: the external 
edge carries the time T^_ 6 , the first graph contains one internal edge with the labeling (0, 0) whereas the second 
graph contains one half-edge with the labeling 2. 
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Appendix A Calculating (r3 5 _3T2} ff 

The Kontsevich model partition function Z-& = e-LtdA*,] satisfies the set of Virasoro conditions: 

L n Z K = 0, n = —1,0,1,..., (A.l) 


where 


L n — (2 n + 3)!! 


d 


dX : 


2 (n+l) 


, (2* + 2n + 1)!!^ d 

+ --TTTi- 


i =0 


(2 i - 1)!! 


n— 1 


+-J^(2n-2i-l)!!(2z + l)!! 


d 2 


i =0 


dT 2 idT 2 ( n -i- 1 ) 


9T 2 (i +n ) 

+ ^11,-iTq /2 + 5n t o 1/16. 


(A.2) 


These conditions for n = — 1 and n = 0 are the corresponding string and dilaton equations which give 
conditions on lower ^-classes: 




S 

( 2 5 -2 + s)^JJr r ^ . 

i =1 ® 


(A.3) 

(A.4) 


For Li we have 


Li = -3- 5 


d 

BTa 


+E(» + 3)(2i + l)T 2i 

i=0 


<) 

dT 2 i +2 


1 d 2 
+ 2 dT 2 


(A.5) 








and gathering the coefficients of the linear term in Tq 9 -q, we obtain 


3 • 5 (t 2 T 3 5 _ 3 ) 9 = (65 - 3)(6c/- 5)(r 3g _ 2 ) 9 + ^(roT 3 9 _ 3 ) 9 _i 

= [3(2 5 - 1)(6 5 - 5) + 2Ag/2](r 3g . 2 ) a = 3[12g 2 - 125 + 5](t 39 _ 2 ) 9 , 

We therefore obtain the following technical lemma 
Lemma A. 9. For the intersection indices, we obtain 

{ T 2 T 3g - 3 )g = ^[12 g{g - 1) + 5](t 39 _ 2 ) s , g > 2. (A. 6 ) 
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